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Abstract 

It is shown that the presence of mixed-culture growth in batch fermentation pro- 
cesses can be very accurately inferred from total biomass data by means of the 
wavelet analysis for singularity detection. This is accomplished by considering sim- 
ple phenomenological models for the mixed growth and the more complicated case 
of mixed growth on a mixture of substrates. The main quantity provided by the 
wavelet analysis is the Holder exponent of the singularity that we determine for 
our illustrative examples. The numerical results point to the possibility that Holder 
exponents can be used to characterize the nature of the mixed-culture growth in 
batch fermentation processes with potential industrial applications. Moreover, the 
analysis of the same data affected by the common additive Gaussian noise still lead 
to the wavelet detection of the singularities although the Holder exponent is no 
longer a useful parameter. 
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1 Introduction 



The growth of microbial species in media containing two or several growth- 
limiting substrates is of great importance in biotechnology and bioengineering. 
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The mixed-culture growth occurs in many industrial processes. A first signifi- 
cant class of such processes is the traditional fermented foods and beverages in 
which either endemic microorganisms or an inoculum with selected microor- 
ganisms are used, see for instance (jH). Some beverages get two or more different 
microorganisms in the inoculum with the purpose to provide a desired flavor. 
Evidence of this influence are presented in the recent paper of (ji]), in which 
the role of different yeast interactions on the wine flavor is discussed. However, 
the phenomenological details and the theory of the time evolution of the fer- 
mentation are as yet poorly understood. We can also mention the interesting 
case of the bioethanol production, in which the substrates used for fermenta- 
tion typically consist of a mixture of glucose and fructose. Bioethanol is the 
product obtained from the metabolism of microbe mixtures feeding with this 
combination of hexoses and pentoses, see e.g., (Q). The last relevant example 
we give is bioremediation, in which gasoline and chemical spills generally yield 
a complex mixture of water-soluble organic compounds. In gasoline spills, for 
instance, the four basic compounds are benzene, toluene, ethylbenzene, and 
xylene. The consumption of this mixture by microorganisms is what is defined 
as the bioremediation process. 

In all the aforementioned cases, the presence of different populations of mi- 
croorganisms and substrates is a key factor in the quality and quantity of the 
final product. Therefore, it is quite useful to detect the presence or lacking 
of process of mixed-culture growth type. Their presence could be used as an 
estimate of the right evolution of the process in its early stage. In addition, a 
rapid and reasonably accurate test is always useful for saving time and helping 
to take quick decisions. It is quite clear then that the biomass concentration 
is one of the most needed quantity that should be measured in fermentation 
monitoring. The most popular method to get the biomass concentration is by 
means of the measurement of the optical density of centrifugalized samples. 
However, this procedure has limited usefulness because it cannot distinguish 
neither the living cells from the dead ones, nor the different types of microor- 
ganisms involved in the process. In some cases it is also possible to correlate 
the total biomass concentration with the values of the redox potential of the 
fermentation. 

Recently, new techniques have emerged to quantify the biomass and distin- 
guish the different microorganisms present in a mixed-culture. Some of then 
based on sophisticated equipment ((0), (@) and (0)) and others resides on 
molecular biology techniques ((0) and (jj)). All these techniques are very 
promising in the study of the dynamics of the mixed-culture growth, although, 
they require expensive or complicated procedures. In this paper, we show that 
it is possible to infer mixed-culture growth of microorganisms from their to- 
tal biomass data, without using such complicated techniques. The alternative 
procedure that we put forth here is based on treating the total biomass data 
by means of the wavelet approach for detection of singularities in the growth 



2 



curves. The idea is to treat the mixed growth curves as more or less regu- 
lar signals that can nevertheless display singularities due to their compound 
structure. In the wavelet literature there exist fundamental papers in which it 
has been shown that the wavelet techniques are very efficient in detecting any 
type of singularities. 

The rest of the paper is organized as follows. In Section 2, we introduce a simple 
dynamics of the mixed-growth type and discuss its basic assumptions. Next, in 
Section 3, the method of the wavelet singularity analysis is briefly presented, 
whereas its application to the mixed type dynamical curves is enclosed in 
Section 4. A conclusion section ends up the paper. An appendix containing 
the standard definitions of Holder exponents of singularities of functions is 
included as well. 



2 A simple mixed-growth model 



The technology of batch processes is well developed and numerous products 
are obtained in this way. Some products such as food, beverages, and pharma- 
ceutical ones require precise tracking of the batch information for safety and 
regulatory purposes. The primary objective of monitoring batch processes is 
to ensure that significant and sustained changes in the quality of the prod- 
uct (caused by disturbances and/or faults) are detected as soon as possible. 
In that sense, the rapid detection of singularities in the output of the batch 
processes offers an interesting solution. The wavelet analysis for singularity 
detection is by now well established but there was no direct application to 
infer mixed-growth in the case of batch biochemical processes. 

In order to achieve this task we will consider here a fermentation process con- 
sisting of a perfectly stirred tank, where no streams are fed into it. In the batch 
fermenter the substrate is converted by biomass into additional biomass and 
products. The general unstructured mass balances for the well-mixed biore- 
actor can be represented by the following equations for the concentrations of 
the cells and substrates: 



dx h 
dt 



dX2i X\i , . 

-w= -if m* 2 ,o 

where X\^ represent the biomass concentrations, Xi,i substrate concentrations 
and Yi is the biomass yield, /Xj (a^.j) is the specific growth rate and i e 1 
represent the i-th species, allowing for the possibility of multiple kinds of 
substrates and microorganisms. The growth rate relates the change in biomass 
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concentrations to the substrate concentrations. Two types of relationships for 
fii (x2 t i) are commonly used: the substrate saturation model (Monod Equation) 
and the substrate inhibition model (Haldane Equation). Both cases will be 
treated here. The substrate inhibited growth can be described by 



where K\ i is the saturation (or Monod) constant, K 2 is the inhibition constant 
and fj, maXl is the maximum specific growth rate. The value of expresses the 
affinity of biomass for substrate. The Monod growth kinetics can be considered 
as a special case of the substrate inhibition kinetics with K 2i = when the 
inhibition term vanishes. For the sake of simplicity, we will consider only two 
species and two substrates. Moreover, we consider that it is possible to measure 
only the total biomass concentration. That means that the output of the 
system (y) will be given by 



where m is the number of species of microorganisms growing in the bioreactor 
(in this work m = 2). We focus on the following four cases: 

I The microorganism and substrate concentrations have the same initial con- 
ditions, but different growth rates, one with a Haldane type and one with a 
Monod type. In addition, quite different values of the Monod constant will 
be taken into account. 
II The microorganism and substrate concentrations have different initial con- 
ditions, but the same growth rates. 

III The microorganism and substrate concentrations have different initial con- 
ditions and different growth rates, one with a Haldane type and one with a 
Monod type. 

IV The microorganism and substrate concentrations have the same initial con- 
ditions and the same growth rates, but with different values of the maximal 
growth rate. 

Table 1 shows the variables and parameter values used to simulated the two 
species growing in the two different substrates, under the four cases under 
consideration. 




m 



y= Y, x ±,i 



i=l 
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Table 1 

The initial conditions and the values of the employed parameters of the mixed- 
growth process model. 



Symbol 


Meaning 


Values 








Units 






Case I 


Case II 


Case III 


Case IV 






Maximal growth rate 


1 


1 


1 


0.9 


11/ h\ 




Saturation parameter 


0.03 


0.03 


0.03 


0.03 


\q/1] 

Li» / "J 




Saturation parameter 


0.5 


0.5 


0.02 


0.5 


\q/l] 

Li/ / "J 


Y\ 

1 


Yield coefficient 


0.5 


0.5 


0.5 


0.5 




x ii 


Initial biomass cone. 


0.1 


0.1 


0.1 


0.25 


\a/l] 

Li// 'J 


•*'2i 


Initial substrate cone. 


10 


10 


10 


10 


\all\ 

lit 1 b \ 




Maximal growth rate 


1 


1 


1 


1 


\llh\ 


K l2 


Saturation parameter 


0.3 


0.03 


0.03 


0.03 


[9/1] 


K 22 


Inhibition parameter 





0.5 





0.5 


[1/9] 


Y 2 


Yield coefficient 


0.5 


0.5 


0.5 


0.5 






Initial biomass cone. 


0.1 


0.2 


0.1 


0.25 


[9/1] 


7 ,o 


Initial substrate cone. 


10 


5 


6 


10 


[9/1] 



3 Mixed cultures on mixtures of substrates 



When microbes are grown in a batch reactor containing a surplus of two sub- 
strates, one of the substrates is generally exhausted before the other, leading to 
the appearance of two successive exponential growth phases. This phenomenon 
could be noticeable at simple view in the biomass signal or unnoticed due to 
its nature or due to additive noise present in the signal. In general, such type 
of phenomenon is known as growth of mixed cultures on mixtures of substrates 
(MCMS). 

The growth of MCMS is a phenomenon of practical and theoretical interest. 
The fundamental understanding of this problem has impact on many practical 
fields such as food processing, production of ethanol from renewable resources, 
bioremediation and microbial ecology, among many others. 

To study the usage of the wavelet approach in the detection of MCMS growth, 
we consider the recent model proposed by Reeves 2004 (Q), which takes into 
account such type of growth. 

Within this section, the index % will denote the species number, and the index 
j will stand for the substrate number. Thus, q denotes the concentration of 
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the zth species, Sj denotes the concentration of the jth substrate, denotes 
the concentration of the lumped system of inducible enzymes catalyzing the 
uptake and peripheral catabolism of Sj by Cj. Here, and Sj are based on the 
volume of the chemostat, and expressed in the units gdw/1 and g/1, respec- 
tively, ejj is based on the dry weight of the biomass, and expressed in the units 
g/gdw. 



y 4J k s 4- 



ij ij l 3 



X 



'J 



13 11 K\; s ■ r, 



ast 



r ast _ 
ij -^ij 

r d . = K d e- 



ds 
~dt 



i-DU- Sj )-r s ljCl -r s 2j c 2 (1) 



dt !<',-<,;*, " %3 V % 13 { ' 



d T rB C{j a. 



where 



K e 



K ij 



VA 



a 



1.1 



Kl 



+ Sj 



Yn rk + Y l2 r\ 2 



Reeves et al (p) comment that a plausible experimental situation is the case 
of Escherichia coli and Pseudomonas aeruginosa, in which, E. coli prefers a 
sugar over an organic acid, and P. aeruginosa prefers the organic acid over 
the sugar. 

In order to have the batch regime in the bioreactor we set parameter D = 
and also employ Reeves' parameters s{ — 1 and s{ = 2. Table 2 shows the 
rest of the parameter values used to simulate the growth of the two species on 
the two different substrates in the MCMS conditions. 
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Table 2 



Parameter values used in the MCMS growth model (J8|) 



V s — 


1000 


V 12 — 


1000 


V s - 
v 21 — 


1000 


V s — 

v 22 — 


1000 


g/gh 


K s n = 


= 0.01 


"l2 ~ 


= 0.01 


K s 2l = 


= 0.01 


K S 2 2 = 


= 0.01 


g/1 


VA = 


0.0025 


V 12 — 


0.0020 


v 2 \ = 


0.0006 


vi 2 = 


0.0036 


g/gdwh 


K e - 


= 0.0017 


■"■12 ~ 


= 0.0032 


"21 - 


= 0.0013 


"22 - 


= 0.0030 


g/gdw 


h.d _ 
K n — 


0.01 


h.d . 
ft 12 _ 


0.01 


h.d . 
ft 21 — 


0.01 


Ud _ 
^22 — 


0.01 


1/h 


k* — 
K U — 


10- 2 V u 




10- 2 Via 


K 21 — 


10- 2 V 2 1 


ft 22 — 


10-2 V 22 


g/gdwh 


Y n = 


0.41 


Y 12 = 


0.24 


Y 2 1 = 


0.35 


^22 = 


0.20 


g/g 



4 Measuring regularity with the wavelet transform 



Let us think of the total biomass of the mixed-growth curves as a signal. 
In general, performing the analysis of a signal means to find the regions of 
its regular and singular behavior. Usually the singularities are very specific 
features for signal characterization. As it has been pointed in the seminal paper 
of (JlJ), the regularity of a signal treated as a function can be characterized 
by Holder exponents. The wavelet transform has been demonstrated to be a 
tool exceptionally well suited for the estimation of Holder exponents (for their 
definitions see the Appendix). 



4.1 The wavelet transform 



Let L 2 (]R) denote the space of all square integrable functions on K. In signal 
processing terminology, L 2 (R) is the space of functions with finite energy. Let 
ip(t) G L 2 (]R) be a fixed function. The function ip(t) is said to be a wavelet if 
and only if its Fourier transform, ip(u>) = J e luJt ip(t)dt, satisfies 



duj < 00. 



\U!\ 



(4) 



The non-divergent relation given by Eq. (4) is called the admissibility condition 
in wavelet theory, see for instance (I12l ) and (114l ). It implies that the wavelet 
must have a zero average on the real line 



ip(t)dt = ^(0) = 0, 



(5) 



and therefore it must be oscillatory. In other words, if) must be a sort of 



wave ( (|l2t Il4j)). Based on if)(t), one defines the functions if) atb as follows 




< 6 (t) = -=4[ , (6) 



where b G M is a translation parameter, while a G M + (a 7^ 0) is a dilation or 
scale parameter. The factor a" 1 / 2 is a normalization constant such that ip a ,b 
has the same energy for all scales a. One notices that the scale parameter a in 
Eq. (6) is a measures of the dilations of the spatial variable (t — b). In the same 
way the factor a -1 / 2 measures the dilations of the values taken by if). Because 
of this, one can decompose a square integrable function f(t) in terms of the 
dilated-translated wavelets ipabif)- We define the wavelet transform (WT) of 
f{t) G L 2 (R) by 



W f (a,b) = (f,i/j a , b )= J f(t)ip a>b (t)dt 

—00 

fm {^) dt, (7) 



1 

7/ 



-00 



where ( , ) is the scalar product in L 2 (R) defined as (f,g) '■— J f(t)g(t)dt, 
and the bar symbol denotes complex conjugation. The WT given by Eq. (7) 
measures the variation of / in a neighborhood of size proportional to a centered 
on point b. In order, to reconstruct / from its wavelet transform (7), one needs 



a reconstruction formula, known as the resolution of the identity ( (11 2c Il4j)). 



00 00 

f(t) = ?r I I W f (a,b)ip a , b (t) 

—00 



da db 
~a^' 



From the above equation we can see why the condition given by Eq. 4 should 
be imposed. One fundamental property that we require in order to analyze 
singular behavior is that if>(t) has enough vanishing moments as argued in the 



works of (|9|) and (113l ). A wavelet is said to have n vanishing moments if and 
only if it satisfies 



00 

J t k ip(t)dx = 0, for k = 0, 1, . . . , n - 1 (9) 

—00 

and 
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J t k i)(t)dt ^ 0, for k > n. (10) 



This means that a wavelet with n vanishing moments is orthogonal to poly- 
nomials up to order n— 1. In fact, the admissibility condition given by Eq. (4) 
requires at least one vanishing moment. So the wavelet transform of f(t) with 
a wavelet ip[t) with n vanishing moments is nothing but a "smoothed version" 
of the n— th derivative of f(t) on various scales. In fact, when someone is in- 
terested to measure the local regularity of a signal this concept is crucial (see 



for instance (Il2t Il4f)). 



4-2 Wavelet singularity analysis 



The local regularity of a function / at a point t is often measured by its Holder 
exponent. The Holder exponent a measures the strength of a singularity at a 
particular point to, where to belongs to the domain of /, see the Appendix. 
It is important to point out that if the singular part of a function / in the 
neighborhood of t is of the type \t — t \ a , then it corresponds to a cusp and in 
this case the singular behavior is fully characterized by its Holder exponent. 
However, there exists functions that involve oscillating singularities which have 



to be described by an additional quantity: an oscillating exponent ( (jlOl ; llll)). 
In such the oscillation has to be analyzed carefully. Such functions can 

not be fully characterized only by the Holder exponent. In this work, we will 
only consider functions whose singularities are not oscillating. 

One classical tool to measure the regularity of a function f(t) is to look at the 
asymptotic decay of its Fourier transform f(u) at infinity. However, the Fourier 
transform is not well adapted to measure the local regularity of functions, 
because it is glob al and provides a description of the overall regularity of 



functions ( (1131 ; Il4f)). Consequently, we need another way to characterize local 



signal regularity. 

In the works (j9; 12; 13; 14) it is shown that the WT provides a way of doing a 
precise analysis of the regularity properties of functions. This is made possible 
by the scale parameter. Due to its ability to focus on singularities in the 
signals, the WT is sometimes referred to as 'mathematical microscope' 



where the wavelet used determines the optics of the microscope 



and its magnification is given by the scale factor a. 

The WT modulus maxima (WTMM) decomposition introduced by (jl3l ) pro- 
vides a local analysis of the singular behavior of signals. In the works of Mallat 



(jl3l ; Il4l ) it has been shown that for cusp singularities the location of the sin- 
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gularity can be detected and the related exponent can be recovered from the 
scaling of the WT along the so-called maxima line (Wtmml for short), which 
is convergent towards the singularity. This is a line where the WT reaches 
local maximum with respect to the position coordinate. Connecting such local 
maxima within the continuous WT 'landscape' gives rise to the entire tree of 
maxima lines. Restricting oneself to the collection of such maxima lines pro- 
vides a particularly useful representation of the entire WT. It incorporates the 
main characteristics of the WT: the ability to reveal the hierarchy of (singular) 
features, including the scaling behavior. 



An other key concept, in addition to vanishing moments, used to characterize 
the regularity of a function in terms of Wtmm is given next. Suppose that 
if) has compact support [-C, C}. The cone of influence of ip at point t is the 
set of points (a, b) in the scale-space plane or domain, such that to is m the 
support of if> a ,b(t)- We will denote the scale-space plane or domain of the WT 
as the (a, 6)-plane or the (a, 6)-domain. Since the support of ip((t — b)/a) is 
[b — Ca, b + Co], the point (a, b) belongs to the cone of influence of to if 



|6-*o|<Ca. (11) 

The function f(t) has a Holder exponent a G (k, k + 1) at to, if and only if 
there exists a constant A > such that at each modulus maxima (a, b) in the 
cone defined by Eq. (11) one has 

\W f {a,b)\ < Aa a+1/2 , a^O, (12) 



(see (Il3l; Il4f) ) . Here it is assumed that the wavelet has at least n > a vanishing 
moments. If f(t) is regular at to or, if the number of vanishing moments is too 
small, i.e., n < a, one obtains for a — > a scaling behavior of the type 



\W f (a,b)\ <Aa n+1/2 . (13) 

The scaling behavior of the Wtmml is given in Eq. (12) and can be rewritten 
as follows 



\og\W f (a,b)\ < \ogA+ (a + -) log a. (14) 

The global Holder regularity at t is thus the maximum slope — | of log |W/(a, b)\ 
as a function of log a along the maxima line converging to to- 
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5 Results and discussion 



In this section, we present the results we obtained using the singularity detec- 
tion procedure described in the previous section. The signal to be analyzed, 
f(t) = y, represents the evolution in time of the total biomass concentration 
for the fermentation process described in Section 2 that includes four different 
cases as specified therein. In all the wavelet-related calculations we employed 
as mother wavelet the first derivative Gaussian tp'(t) = d/dt(e~ t / 2 ) having 
only one vanishing moment. The final goal is always to calculate the Holder 
exponent of the singularities for such processes because it is a direct measure 
of the irregularity of a signal (function) at the singular point to, in the sense 
that higher values of it correspond to more regular functions than the lower 
values. 

Figure 1 a,b,c shows the performance of the wavelet singularity analysis as 
applied to Case I (same initial conditions but different kinetic rates). We 
obtain a Holder exponent of quite high value. 

The following figure shows the performance of the scheme applied to Case II 
(same growth rates but different initial conditions). In this case, the Holder 
exponent of the mixed growth singularity is lower than in Case I. 

Similarly to the previous cases, Fig. (3) presents the graphical results for 
Case III (different initial conditions and different growth rates). Although, 
the singularity looks very mild in the time evolution of the total biomass 
concentration the wavelet analysis is able to detect it with high precision. 

Finally, Case IV (same initial conditions, same growth rates but with different 
values of their maximal growth rates) is graphically analyzed in Fig. (4). For 
this case we obtained the lowest Holder exponent. 

Although the latter two cases seem to correspond to almost overlapping of the 
WTMML pointing to bifurcation phenomena we are still not at the threshold 
of a completely different behavior of the log plots generated by bifurcations. 
This could be explained by the fact that the strength of the first singularity 
is bigger with respect to the second one. 

5. 1 Wavelet analysis for the MCMS case 

The MCMS case is the most interesting case that we discuss here because we 
will show that it is possible to infer in a very accurate manner by means of 
WT the moment in which the microorganisms switch their carbon source. In 
order to understand the detailed dynamics of this combined growth, we first 
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apply separately the WT approach to the two biomass signals y = c\ (Fig. 5) 
and y = c 2 (Fig. 6) and then to the total signal y = c\ + c 2 (Fig. 7). 

It is worth noting that the Holder exponent is bigger than one, a quite inter- 
esting feature which means that the singularity lies in the second derivative 
of the biomass signal. This result gives further opportunities to character- 
ize the nature of the singularity because it suggests that the type of growth 
can be inferred from the order of the derivative in which the singularity oc- 
curs directly given by the value of the Holder exponent. The latter fact is a 
great simplification with respect to the analytical search of the singularities 
which implies obtaining the analytical solution of the given dynamical growth 
model. Moreover, even in such fortunate cases, the analytical solutions could 
be subject to fixed parameter values of the model. On the other hand, the 
WT numerical approach allows the singularity analysis even in the case of 
time- varying parameters. 



5.2 Wavelet analysis for the noisy data case 



It is well known that in some cases the amplitude of the Gaussian noise af- 
fecting the on-line signals can be an important annoying factor. Therefore, a 
good analysis should be robust in such cases. Thus, we provide here the WT 
analysis for the MCMS biomass signals in the presence of white noise, that 
is, in the next figures (8-10) we consider signals of the form yi = Ci + e(t) and 
their sum, where e(t) stands for the functional form of the noise. 

We notice from the corresponding plots that the noisy data do not allow to 
obtain global Holder exponents in a straightforward manner, which is a result 
already reported in the wavelet literature ((14) and (flil)). On the other hand, 



the singularity detection is robust with respect to the noise for reasonable 
levels of its amplitude. In addition, Figure 9 gives us the hint that when the 
cones of influence produced by the Gaussian noise enter the scales of the 
singularity the cone of the latter becomes undistinguishable from those of 
the noise. This remark could be used as a sort of resolution criterium of the 
WT method in the presence of noise. Therefore, one can determine a critical 
amplitude of the noise for which the WT approach looses its applicability. 



6 Concluding remarks 



We showed here explicitly how the wavelet singularity analysis can be applied 
to infer mixed growth behavior of fermentation processes using only total 
biomass data. We prove that the wavelet analysis is very accurate for all the 
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cases we considered. A very interesting feature of our research is that the 
Holder exponent is sensitive to the type of the mixed-growth phenomenon, 
more specifically depends on the parameters of the growth processes and on 
their initial conditions. The MCMS case points to the remarkable technological 
possibility of detecting the change of the substrate uptake since the singularity 
appears in the second derivatives of the biomass signal. This can lead one to 
think of the possibility to infer substrate contaminations based only in the 
analysis of the biomass data. In addition, our results for the noisy data clearly 
hint to the fact that the wavelet singularity analysis maintains its attractive 
features even in these more difficult but realistic case. We hope that in future 
works we could find out the mathematical relationships implied by this possible 
correlation. It might allow the usage of the Holder exponent as an identification 
criterium of the more specific nature of mixed-growth processes. 
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Appendix 



A function / : R — > R is said to be Holder continuous of exponent a (0 < a < 
1) if, for each bounded interval (c, d) C R, we can find a positive constant A 
such that 



i/(*)-/(t )i <K\t-t \ a (15) 

for all t, to G (c, d). 



The space of Holder continuous functions is denoted C a . A function is said 
to be C n+a if it is in C n and its nth derivative is Holder continuous with 
exponent a. Thus, if we consider the Holder exponent n < a < n + 1, with 
nGN, the function can be differentiated n times, but the (n + l)th derivative 
does not exist. Therefore, a function with a Holder exponent n < a < n + 1 
is said to be singular in the nth derivative. Keeping this in mind, let us give 



the following definition of the Holder regularity of a function (Il2t Il3l : 



Let nGN and n<a<n + l.A function f(t) has a local Holder exponent 
a at to if and only if there exist a constant A > 0, and a polynomial P n (t) 
of order n, such that 

Vt G R, |/(f) - P n (t - t ) I < K\t - t \ a (16) 

The function f(t) has a global Holder exponent a on the interval (c, d) if 
and only if there is a constant A and a polynomial of order n, P n (t), such 
that equation (16) is satisfied for all t G (c, cf). 

The Holder regularity of /(£) at t is the supremum of the a such that /(t) 
is Holder a at to- 

The nth derivative of a function f(t) is singular at to if f(t) has a local 
Holder exponent a at t with n < a < n + 1. 



A function /(t) that is continuously differentiable at a given point has a Holder 
exponent not less than 1 at this point. If a G (n,n + 1) in (16) then f(t) is n 
times but not (n + 1) times differentiable at the point to, and the polynomial 
P n (t) corresponds to the first [n + 1) terms of the Taylor series of f(t) around 
t = t . For example, if n = 0, we have Po(t — to) = /(to)- 
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Figure Captions 



Fig. 1 

a) The time evolution of the total biomass concentration signal for Case I. 

b) The wavelet cones of influence corresponding to this case showing a very 
accurate identification of the two singularity points presented in the signal, of 
which the first one allows to infer the presence of the mixed growth feature of 
the fermentation process whereas the second one is associated with the end of 
the fermentation batch cycle. 

c) From the slope in the double logarithmic plot, the Holder coefficient of the 
mixed growth singularity is calculated as a = 0.95. 

Fig. 2 

a) The time evolution of the total biomass concentration signal for Case II. 

b) The wavelet cones of influence corresponding to this case again showing 
the accurate identification of the two singularity points, of the same type, 
respectively, as in Fig. 1. 

c) The Holder coefficient of the mixed growth singularity is now a = 0.88. 
Fig. 3 

Same caption comments as in the previous figures but for Case III. The Holder 
coefficient of the mixed growth singularity is now a = 0.92. 

Fig. 4 

Same caption comments as in the previous figures but for Case IV. The value 
of the he Holder coefficient for the mixed growth singularity is a = 0.84. 

Fig. 5 

Same caption comments as in the previous figures but for the MCMS biomass 
signal y — C\. The Holder coefficient of the mixed growth singularity is now 
a = 1.89. 

Fig. 6 

Same caption as in the previous figures but for the MCMS biomass signal 
y = c 2 . The Holder coefficient of the mixed growth singularity is now a = 1.87. 
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Fig. 7 



Same caption comments as in the previous figures but for the total MCMS 
signal. The Holder coefficient of the mixed growth singularity is now a = 1.88. 

Fig. 8 

MCMS data corresponding to c\ with a small amplitude Gaussian noise added. 
From the bottom plot c) one can see that because the curve is not a straight 
line one cannot get a global Holder coefficient from its slope. 

Fig. 9 

MCMS data corresponding to c 2 with a small amplitude Gaussian noise added. 
The Holder coefficient is not a useful concept in this case. 

Fig. 10 

MCMS data corresponding to the sum c\ + c 2 with the same Gaussian noise 
added. The concept of global Holder coefficient is again not useful. 
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This figure "Casl.png" is available in "png" format from: 



http://arXiv.org/ps/physics/05 121 86v2 



This figure "mcmsl.png" is available in "png" format from: 



http://arXiv.org/ps/physics/05 121 86v2 



This figure "Cas2.png" is available in "png" format from: 



http://arXiv.org/ps/physics/05 121 86v2 



This figure "mcms2.png" is available in "png" format from: 



http://arXiv.org/ps/physics/05 121 86v2 



This figure "Cas3.png" is available in "png" format from: 



http://arXiv.org/ps/physics/05 121 86v2 



This figure "Cas6.png" is available in "png" format from: 



http://arXiv.org/ps/physics/05 121 86v2 



This figure "mcmsT.png" is available in "png" format from: 
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This figure "mcmsrl.png" is available in "png" format from: 
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This figure "mcmsr2.png" is available in "png" format from: 
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This figure "sumar.png" is available in "png" format from: 
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